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Abstract 

The  kinetics  of  CO  preferential  oxidation  (PROX)  was  studied  to  evaluate  various  rate  expressions  and  to  simulate  the  performance 
the  CO  oxidation  step  of  a  methanol  fuel  processor  for  fuel  cell  applications.  The  reaction  was  carried  out  in  a  micro  reactor  testing  unit 
using  a  commercial  Engelhard  Selectoxo  (Pt-Fe/y-alumina)  catalyst  and  three  self-prepared  catalysts.  Temperature  was  varied  between 
100  and  300  °C,  and  a  of  range  feed  rates  and  compositions  were  tested.  A  reaction  model  in  which  three  reactions  (CO  oxidation,  H2 
oxidation  and  the  water  gas  shift  reaction)  occur  simultaneously  was  chosen  to  predict  the  reactor  performance.  Using  non-linear  least 
squares,  empirical  power-law  type  rate  expressions  were  found  to  fit  the  experimental  data.  It  was  critical  to  include  all  three  reactions  to 
determine  good  fitting  results.  In  particular,  the  reverse  water  gas  shift  reaction  had  an  important  role  when  fitting  the  experimental  data 
precisely  and  explained  the  selectivity  decrease  at  higher  reaction  temperatures.  Using  this  three  reaction  model,  several  simulation  studies 
for  a  commercial  PROX  reactor  were  performed.  In  these  simulations,  the  effect  of  O2/CO  ratio,  the  effect  of  water  addition,  and  various 
non-isothermal  modes  of  operation  were  evaluated.  The  results  of  the  simulation  were  compared  with  corresponding  experimental  data 
and  shows  good  agreement. 
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1.  Introduction 

Removing  carbon  monoxide  in  a  hydrogen  rich  stream 
is  a  critical  issue  and  an  unavoidable  problem  when  hydro¬ 
carbons  are  used  as  the  hydrogen  source  for  fuel  cells.  Re¬ 
gardless  of  the  reformer,  small  amounts  of  CO  exist  after 
hydrocarbon  reforming  and  even  after  the  water  gas  shift 
reaction.  These  small  amounts  of  CO,  typically  less  than 
1  mol%,  must  be  removed  to  prevent  poisoning  of  the  fuel 
cell  electrodes.  Among  the  various  methods  to  remove  CO 
selectively,  catalytic  oxidation  is  considered  as  one  of  the 
most  plausible  and  economical  options.  Due  to  the  recent 
growth  in  research  on  fuel  cells  and  fuel  processing,  a  large 
number  of  studies  on  CO  selective  or  preferential  oxidation 
(PROX)  have  been  published. 

The  catalysts  for  CO  selective  oxidation  can  be  classified 
into  three  categories  as  shown  in  Fig.  1,  which  summarizes 
30  recent  publications  [1-30].  The  most  commonly  used  for¬ 
mulation  is  platinum  or  other  precious  metals  on  alumina, 
and  at  temperatures  around  200  °C.  Also,  Au  based  cata¬ 
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lysts  show  good  performance  at  lower  temperatures  around 
100  °C  which  is  close  to  PEM  fuel  cell  operating  temper¬ 
atures.  Lastly,  several  common  transition  metals  have  been 
investigated  to  find  a  more  economical  CO  selective  oxida¬ 
tion  catalyst.  Most  recent  papers  are  focus  on  catalyst  for¬ 
mulation,  characterization,  and  basic  performance  such  as 
activity  and  selectivity  of  CO.  Few  papers  have  investigated 
the  kinetics  and  rate  expressions  of  the  reactions  involved. 

In  the  CO  selective  oxidation  reaction  system,  the  follow¬ 
ing  two  oxidation  reactions  (1)  and  (2)  occur. 

CO  +  \  02  -*  CO2,  A H°29%  =  -283  kJ/mol  (1) 

H2  +  \  02  -*  H20,  A H°29S  =  -242  kJ/mol  (2) 

For  precious  metal  catalysts,  several  different  kinetic  rate 
expressions  have  been  reported.  Amphlett  et  al.  [4]  devel¬ 
oped  a  simple  first-order  rate  expression  for  CO  oxidation 
on  a  platinum/alumina  catalyst  as  follows: 


f'co  =  /ccocco(mol s  1  kg 


where  kco  =  0.026  exp 


(m3  s  1  kg  l) 
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Co,  Fe,  Cu, 
17.0% 


Au  based, 
13.2% 


Others  (Mn, 
Cr,  Mg,  Zn), 
9.4% 


Pt  based, 
35.8% 


Ru,  Rh,  Pd, 
24.5% 


Fig.  1.  Typical  formulation  of  CO  selective  oxidation  catalysts  in  the  30 
recently  published  papers. 


Kahlich  et  al.  [7]  derived  a  different  kinetic  expression  by 
introducing  a  process  parameter  A  which  is  defined  as  the 
concentration  ratio  of  oxygen  to  CO  (Eq.  (4)).  According 
to  this  kinetic  expression,  the  reaction  rate  of  CO  can  be 
determined  with  reaction  orders  of  —0.42  for  CO  and  +0.82 
for  O2  and  an  activation  energy  of  71  kJ  mol-1. 


j  p0.42  0.82 

r CO  =  M  Fco  A 


where  A  = 


2Co9  2/\> 


Cco  Pco 


(4) 


kcK^PcoP1^ 
1  +  ^P'co 


In  spite  of  the  many  papers  that  give  suggested  rate  expres¬ 
sions  as  described  above,  there  are  very  few  papers  that 
consider  H2  oxidation  simultaneously  with  CO  oxidation 
reaction  in  PROX  related  publications. 

Rate  expressions  for  hydrogen-oxidation  reaction  over  the 
platinum  catalyst  have  been  found  separately  in  several  pub¬ 
lications  [31-33].  The  rate  expressions  of  power  law  type 
empirical  expressions  are  summarized  in  the  Table  lb.  As 
shown  in  these  rate  expressions,  the  H2  oxidation  reaction 
rate  mainly  depends  on  the  partial  pressure  of  oxygen  and 
not  the  partial  pressure  of  hydrogen. 

The  feed  stream  of  a  PROX  reactor  in  a  methanol  fuel  re¬ 
former  system  is  composed  of  H2,  O2,  CO,  CO2  and  H2O. 
Therefore  in  addition  to  the  reactions  of  CO  and  H2  oxida¬ 
tion,  it  is  also  necessary  to  consider  the  reverse  water  gas 
shift  reaction  in  the  kinetic  model  of  the  PROX  reactor. 


Also,  Kim  and  Lim  [20]  used  a  similar  rate  expression  as 
follows: 


— ^co,tof  =  1.4  x  108exp 


-78 

~RT 


p—0.51  p0.76 

rco  ro2 


In  Table  1,  the  kinetic  parameters  of  several  empirical  power 
law  rate  expressions  are  compared. 

For  another  approach,  Venderbosch  et  al.  [9]  (Eq.  (6)) 
found  the  parameters  of  the  Langmuir-Hinshelwood  (LH) 
expression  using  a  downhill-simplex  minimization  method 
and  experimental  data. 


r(c ,  7)  = 


4’  exp(-Ea/RT)cCoco2 

(1  +  K®  exp(A//ads/7?T)cCO)2 


Finally,  Sedmak  et  al.  [27]  evaluated  two  LH  type  rate  ex¬ 
pressions  for  the  Cu-CeO  based  catalyst  shown  in  Eqs.  (7) 
and  (8). 

kcoko2  Pco  Pq? 

rco  =  - - - -  (7) 

0.5Aco^co  +  ko2  P0l 


H20  +  CO  C02  +  H2,  A 7/298  =  -41.1  kJ/mol 

(9) 

The  equilibrium  constant  of  the  water  gas  shift  reaction  is 
about  225  at  200  °C  suggesting  the  forward  reaction  should 
be  dominant  in  this  temperature  region  if  the  gas  does  not 
have  large  amounts  of  H2and  CO2.  However,  in  most  PROX 
reactor  feed  gases,  the  concentration  of  H2  and  CO2  are 
approximately  70  and  20%,  respectively,  which  are  high 
enough  to  influence  the  reverse  WGS  reaction.  Therefore,  to 
accurately  predict  the  concentration  of  all  gas  components, 
all  three  reactions  (CO  oxidation,  H2  oxidation  and  (reverse) 
water  gas  shift  reaction)  must  be  considered  simultaneously. 

This  three  reaction  system  and  the  exact  rate  expressions 
for  each  reaction  are  important  components  in  the  optimiza¬ 
tion  and  control  of  commercial  fuel  reformers.  The  goals 
of  this  study  are  to  determine  the  kinetic  models,  to  find 
the  exact  rate  parameters,  and  to  simulate  the  CO  selective 
oxidation  system.  Through  this  approach,  we  hope  to  gain 


Table  1 


Empirical  expressions  for  CO  and  1+  oxidation  in  the  literature 


Catalyst 

a 

p 

ln(*o) 

E  (kJmol  l) 

Reference 

(a)  For  CO  oxidation  rate: 

r  =  koexp(-E/RT)K  P^0 

Pt/y-AbOa 

0.8 

-0.4 

72-76 

[7] 

1 

-1.5 

— 

55 

Muraki  [7] 

0.7 

-0.1 

67-71 

Sarkany  [7] 

Pt/Si02 

1 

—  1  to  0 

56-80 

Au/a-Fe203 

1 

1 

31 

[11] 

Cuo.iCeo.902-y 

Y 

11.88 

57.2 

[27] 

(b)  For  H2  oxidation  rate: 

r  =  k0exp(-E/RT)P^P^2 

Pt/y-A1203 

0 

0.8 

21.9 

[31] 

Pt/y-A1203 

0 

1 

— 

20.9-58.8 

[32] 

Pt(l,l,l) 

— 

0.59 

— 

14.3-47.0 

[33] 
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important  beneficial  insights  for  the  design  and  operation  of 
the  PROX  reactor  in  a  fuel  cell  system. 


2.  Experimental 

2.7.  Catalysts  and  reactors 

For  the  experiments  in  this  study,  a  commercial  catalyst 
and  several  in  house  prepared  catalysts  have  been  tested. 
Reaction  conditions  are  chosen  similar  to  the  reactor  con¬ 
ditions  which  would  be  found  in  a  fuel  processor  system, 
such  as:  feed  gas  flow  rate  and  composition,  space  velocity, 
pressure,  and  reactor  temperature.  The  commercial  catalyst 
used  was  the  Selectoxo  catalyst  manufactured  by  Engelhard 
(Pt-Fe/Al203).  The  Selectoxo  catalyst  is  a  3  mm  x  4  mm 
sized  pellet  catalyst  composed  of  0.5%  Pt  and  0.02%  Fe  on  a 
y-alumina  support.  The  catalysts  pellets  were  ground  lightly 
to  a  particle  diameter  of  100-250  jam  to  eliminate  internal 
diffusion  resistance  but  allow  good  gas  distribution. 

Three  different  in-house  catalysts,  Pt,  Rh,  and  Pt-Rh 
over  the  y-alumina  were  prepared.  The  platinum  cata¬ 
lyst  was  made  from  an  aqueous  solution  of  hydrogen 
hexachloro-platinate(IV)  hydrate  (Aldrich  Chemical  Co.), 
and  the  rhodium  catalyst  was  made  from  rhodium(III)  chlo¬ 
ride  hydrate  (Aldrich  Chemical  Co.)  both  by  the  conven¬ 
tional  impregnation  method.  The  amounts  of  metal  loading 
after  the  impregnation  were  respectively  0.55%  Pt,  0.74% 
Rh,  and  0.55%  Pt-0.07%  Rh.  After  the  impregnation,  all 
three  catalysts  were  dried  at  120  °C  for  12  h  and  calcined 
at  450  °C  for  3h.  Because  the  CO  and  H2  oxidation  re¬ 
actions  are  highly  exothermic,  0.5  g  of  each  catalyst  was 
diluted  with  5.0  g  of  inert  alumina  to  maintain  isothermal 
conditions. 

All  reaction  tests  were  performed  in  a  standard  catalyst 
test  unit.  A  stainless  steel  (or  glass)  tubular  reactor,  1/2  in. 
in  diameter  and  12  in.  long  was  used  for  all  reaction  tests. 
To  ensure  isothermal  conditions  along  the  bed  length,  a  split 
tubular  furnace  was  used  and  the  temperature  of  the  catalyst 
bed  was  measured  directly  by  a  1/16  in.  J-type  thermocouple. 

2.2.  Conditions  and  product  analysis 

To  simulate  the  conditions  exiting  the  water  gas  shift  re¬ 
actor  in  a  methanol  reforming  unit,  the  feed  gas  stream  con¬ 
tained  62-72%  H2,  0. 5-5.0%  02,  2-17%  N2,  0.5-3.0%  CO, 
and  20-24%  CO2  (dry  basis).  In  cases  when  water  is  added, 
the  injection  rates  were  controlled  precisely  by  a  syringe 
pump,  74900  Series  (Cole  Palmer),  from  0.5  to  8mlh_1. 
The  reaction  tests  were  performed  at  temperatures  between 
100  and  300  °C.  With  these  catalyst  loadings  and  feed  com¬ 
positions,  the  GHSV  at  reaction  temperature  was  controlled 
between  1000  and  20,000  h  .  All  the  reaction  runs  were 
performed  under  atmospheric  pressure. 

The  effluent  of  the  reactor  was  connected  directly  to  a 
CARLE  Series  S  gas  chromatograph  containing  a  thermal 


Table  2 


Feed  composition  and  catalysts  in  the  experiment 


Run  No. 

Feed  mixture  (mol%) 

Catalysts 

Water 

CO 

O2 

1-3 

~3% 

3%,  5%,  1.5% 

Pt-Fe  (Selectoxo) 

— 

4 

~1% 

1% 

Pt-Fe 

— 

5-6 

~0.6% 

1.2%,  0.6% 

Pt-Fe 

— 

7-9 

~1% 

1% 

Pt  0.55,  Rh,  Pt-Rh 

— 

10-11 

~1% 

1% 

Pt  0.19%,  Pt  1.6% 

— 

12-13 

~1% 

1% 

Pt-Fe 

15%,  7% 

conductivity  detector  using  helium  as  the  carrier  gas.  The 
column,  a  Supelco  Carboxen  1000  (60-80  mesh,  15  ft  x 
1/8  in.),  was  used  to  analyze  the  light  gas  components.  The 
GC  chamber  was  maintained  at  a  constant  temperature  of 
50  °C  to  separate  N2  and  O2  properly.  Five  components  H2, 
O2,  N2,  CO  and  CO2  were  measured  during  each  test  run. 
Material  balances  on  carbon  were  calculated  to  verify  mea¬ 
surement  accuracy,  and  for  all  runs  reported  here  were  within 
3%  of  closure. 


2.3.  Activity  and  selectivity  measurements 


CO  activity  and  selectivity  (moles  of  CO  oxidized  divided 
by  two  times  the  moles  of  O2  consumed)  were  measured  at 
various  levels  of  CO,  O2,  and  water.  The  catalysts  and  feed 
streams  of  14  runs  are  summarized  in  Table  2.  The  oxygen 
level  (air  addition  rate)  is  an  important  factor  in  PROX  op¬ 
eration.  Theoretically  only  0.5  mol  of  oxygen  is  needed  to 
remove  one  mole  of  CO.  However,  most  of  PROX  reactors 
use  excess  air  to  some  extent  due  to  H2  oxidation  reaction. 
Because  CO  concentrations  entering  a  commercial  PROX 
reactor  are  usually  less  than  1%,  we  varied  the  CO  concen¬ 
tration  from  0.5  to  3.0%  and  varied  the  O2  concentration 
from  0.5  to  6.0%. 

In  each  run,  five  to  eight  different  temperatures  (under 
300  °C)  were  selected.  At  temperatures  under  150  °C,  it  was 
difficult  to  get  repeatable  data  because  of  low  conversion  and 
temperature  fluctuations  between  the  catalyst  surface  and  the 
reactor.  Therefore,  most  of  the  data  in  this  low  temperature 
range  were  obtained  by  averaging  12-24h  of  data  at  the 
targeted  temperature.  CO  conversion  was  calculated  directly 
from  the  concentration  change  of  CO.  The  selectivity  of  CO 
in  the  PROX  reactor  was  defined  as  follows. 


Selectivity  (%)  = 


x  100 


3.  Experimental  results 

3.1.  CO  activity  and  selectivity  over  a  commercial  catalyst 

Fig.  2(a)  and  (b)  show  a  plots  of  CO  conversion  and  se¬ 
lectivity  at  various  O2/CO  feed  ratios  as  a  function  of  re- 
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Fig.  2.  Activity  and  selectivity  for  CO  selective  oxidation  with  various 
O2/CO  ratio  as  a  function  of  reaction  temperature:  0.5  g  of  Selectoxo 
(Pt-Fe/y-alumina)  catalyst,  flow  rate  =  167  seem,  H2  =  64-75  mol%, 
O2/CO  in  mol%,  1  atm,  no  water  addition. 

action  temperature  for  the  Selectoxo  catalyst.  As  shown  in 
this  figure,  the  conversion  of  CO  increases  at  high  O2/CO 
ratios  and  low  CO  concentrations.  When  the  O2/CO  ratio 
was  2.0,  the  maximum  CO  conversion  was  more  than  95% 
even  under  200  °C  (see  R2,  R5).  In  the  three  cases,  where 
the  O2/CO  mole  ratio  is  near  1 .0,  lower  CO  concentrations 
gave  higher  CO  conversions  (see  R4,  R6).  However,  selec¬ 
tivity  was  low  when  the  O2/CO  ratio  was  high  (see  R2,  R5). 
When  the  O2/CO  ratio  is  around  one  or  less,  the  selectiv¬ 
ity  of  the  catalyst  has  a  maximum  value  between  200  and 
250  °C  (see  Rl,  R3). 

3.2.  Other  precious  metal  PROX  catalysts 

The  performance  of  three  in-house  catalysts  is  compared 
with  the  commercial  Pt-Fe  catalyst  in  Fig.  3.  Rhodium  is 
a  common  catalyst  used  for  oxidation  reactions  and  one  of 
the  candidates  for  the  CO  PROX  reaction.  Although  Rh  is 
used  for  automotive  catalysts,  there  are  very  few  papers  ad¬ 
dressing  CO  selective  oxidation  for  Pt-Rh  bimetallic  cata¬ 
lysts.  Cai  et  al.  [34]  reported  the  synergist  effects  of  Pt-Rh 
alloy  particle  at  certain  conditions  for  CO  oxidation.  How¬ 
ever,  that  work  was  when  hydrogen  was  not  a  feed  compo¬ 
nent.  As  shown  in  Fig.  3,  the  rhodium  catalyst  has  a  lower 
activity  than  any  of  the  Pt  catalysts.  The  bimetallic  catalyst 


Fig.  3.  Activity  and  selectivity  for  CO  selective  oxidation  with  various  cat¬ 
alyst  formulations  as  a  function  of  reaction  temperature:  0.5  g  of  catalysts 
(0.55%  Pt,  0.78%  Rh,  0.55%  Pt-0.07%  Rh,  and  0.5%  Pt-0.02%  Fe),  flow 
rate  —  167  seem,  H2  —  72.0  mol%,  O2  -  1.0mol%,  CO  —  1.0mol%, 
CO2  =  22mol%,  1  atm,  no  water  addition. 

(Pt  0.55  wt.%-Rh  0.07  wt.%)  shows  no  synergetic  effect  and 
gives  almost  the  same  results  as  the  Pt  (0.55  wt.%)  catalyst. 
The  Pt-Fe  catalyst  has  the  highest  CO  conversion  and  the 
rhodium  catalyst  shows  the  lowest  selectivity  values.  The 
only  surprising  result  in  Fig.  3  is  that  the  platinum  cata¬ 
lyst  shows  a  slightly  better  selectivity  at  temperatures  above 
250  °C. 

Fig.  4  shows  the  activity  of  the  in  house  platinum  catalyst 
with  varying  Pt  loading  (1.6,  0.55  and  0.19%  Pt).  As  ex¬ 
pected,  the  activity  directly  increases  with  platinum  amount 
over  this  range  of  Pt  loadings.  However,  the  high  loading 
platinum  catalyst  (1.6%)  shows  the  poorest  selectivity  at 
higher  temperatures. 

3.3.  Water  addition  effect 

Water  content  is  one  of  the  most  critical  factors  in  the 
PROX  reaction  and  is  very  important  for  proper  operation  of 
the  fuel  cell  stack.  Contradictory  viewpoints  still  exist  as  to 
whether  water  acts  positively  for  CO  oxidation  or  not.  Many 
research  groups  [8,21,22]  have  reported  that  water  enhances 
the  CO  oxidation  reaction  for  a  certain  range  of  tempera¬ 
tures.  On  the  other  hand,  Korotkikh  and  Farrauto  reported 
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Fig.  4.  Activity  and  selectivity  for  CO  selective  oxidation  with  various 
platinum  loadings  as  a  function  of  reaction  temperature:  0.5  g  of  catalysts, 
platinum  loading  in  weight  %,  flow  rate  =  167  seem,  H2  =  72.0  mol%, 
O2  —  1.0mol%,  CO  —  1.0mol%,  CO2  -  22mol%,  1  atm,  no  water 
addition. 


Fig.  5.  Activity  and  selectivity  for  CO  selective  oxidation  with  addi¬ 
tion  of  water  as  a  function  of  reaction  temperature:  0.5  g  of  Selec- 
toxo  (Pt-Fe/y-alumina)  catalyst,  flow  rate  =  167  seem,  H2  =  72.0  mol%, 
O2  =  1.0mol%,  CO  —  1.0mol%,  CO2  =  22mol%,  1  atm. 


that  the  addition  of  small  amounts  of  water  (3mol%)  de¬ 
creased  the  CO  conversion.  Grisel  and  Nieuwenhuys  [19] 
and  Avgouropoulos  et  al.  [24]  also  reported  that  H2O  had  a 
detrimental  effect  on  the  CO  oxidation  activity  of  Au  based 
catalyst. 

In  this  study,  as  shown  in  Fig.  5,  the  activity  and  selectivity 
increases  with  water  addition  up  to  the  temperature  of  220  °C 
for  water  addition  of  15%  or  7%  (mol%).  Above  220  °C, 
large  amounts  of  water  had  a  positive  effect  on  selectivity 
and  activity  while  smaller  amounts  of  water  were  detrimental 
compared  to  the  dry  case.  The  reasons  for  these  water  effects 
are  discussed  in  the  following  section. 

4.  Kinetics  and  simulations 

4.1.  Kinetics  and  rate  expressions  for  the  PROX 
reactions 

As  described  previously,  this  study  is  distinctive  in  that  it 
addresses  all  three  reactions  in  the  CO  oxidation  system:  H2 
oxidation,  water  gas  shift  reaction,  and  CO  oxidation.  For 
the  water  gas  shift  reaction  (3),  the  kinetics  is  well  known 


and  it  is  reliable  to  write  a  rate  expression  as  follows  [35]: 

7  f  n  „  ^co2^h2\ 

-r3  =  £3  (  Fco^h2o - — — J  , 

/  4577.8  \ 

where  Kp  =  exp  (  — — - 4.33  j  (10) 

The  H2  oxidation  reaction  (2)  in  a  hydrogen  rich  environ¬ 
ment  can  be  modeled  using  an  empirical  rate  expression 
similar  to  others  [31-33]. 

-r2=k2Pf  (11) 

And  an  empirical  expression  can  be  adequately  used  for  the 
rate  of  CO  oxidation  (1). 

-n  (12) 

The  constants  in  each  rate  expression  are  assumed  to  be  Ar¬ 
rhenius  functions  of  temperature  (A0  exp(— E/RT)),  which 
gives  nine  parameters  (three  frequency  factors,  three  acti¬ 
vation  energies,  and  two  power  constants)  to  be  found  for 
the  complete  model.  All  forty  data  sets  shown  in  previous 
graphs  were  used  to  find  these  nine  parameters. 

By  using  data  at  constant  temperature,  the  values  of  a,  ft, 
y,  k\ ,  /c2,  and  k3  were  determined  by  minimizing  the  sum  of 
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the  square  of  the  difference  between  predicted  and  observed 
CO  conversion  (Eq.  (7)). 


Minimize 


N 


y^(Xeco,/  -  Xcco,0' 


i= 1 


(13) 


where  i  =  experiment  number  (N  =  28  for  CO  selec¬ 
tive  oxidation),  Xeco ,i  =  experimental  CO  conversion  at 
the  reactor  outlet  (mol/h),  and  Xeco,/  =  calculated  CO  con¬ 
version  at  the  reactor  outlet  (mol /h) . 

This  procedure  found  that  the  best  power  law  exponents 
were  0.5  for  O2  and  —0.1  for  CO.  Fixing  the  exponents  at 
these  values,  all  of  data  (40  sets)  were  fit  simultaneously  by 
varying  the  activation  energies  and  pre-exponential  factors 
to  again  minimize  the  sum  of  squares  error  (Eq.  (7)).  To 
obtain  the  exiting  concentrations  for  each  experiment,  the 
MATLAB  subroutine  function  ODE23  was  used  for  numer¬ 
ical  integration.  To  minimize  Eq.  (7),  the  MATLAB  subrou¬ 
tine  function  LSQNONLIN  was  used  for  non-linear  least 
squares  optimization. 

The  empirical  rate  expressions  derived  from  this  numeri¬ 
cal  analysis  are  as  follows: 

-n  =  3.528  x  102  exp  (■ - — J  P~ g-1  (14) 

-r2  =  2.053  x  lOexp  ( J  P^5  (15) 


— r3  =  4.402  x  103  exp 

X  (/-con.,0  -  (16) 

The  quality  of  the  fit  is  shown  in  Fig.  6  by  comparing  the 
observed  and  calculated  CO  conversion  and  selectivity  for 
all  the  experiments.  As  shown  in  the  figure,  the  calculated 
values  show  good  agreement  with  the  experimental  values. 
The  R 2  values  for  CO  conversion  was  0.92  for  all  data  sets. 
Although  selectivity  was  not  specifically  fit,  the  calculated 
values  from  the  simulation  are  located  within  ±10%  errors 
range  for  most  experimental  data. 

As  shown  in  Table  1 ,  literature  values  of  the  power  con¬ 
stants  for  CO  oxidation  are  0. 7-1.0  for  oxygen,  and  —1.5 
to  0  for  CO,  and  the  activation  energy  for  CO  oxidation  is 
between  55  and  80kJmol_1.  The  rate  parameters  from  this 
study  have  a  slightly  lower  value  for  02  order  and  activation 
energy.  To  address  this  concern,  the  experimental  data  was 
fit  using  only  one  reaction  (CO  oxidation).  In  that  approach, 
the  resulting  best- fit  power  constant  for  oxygen  was  0.76, 
which  is  similar  to  the  values  in  the  literature,  and  a  higher 
activation  energy  of  46  kJ  mol-1  was  found.  However,  with¬ 
out  using  all  three  reactions,  the  simulated  values  for  CO 
conversion  did  not  compare  well  with  the  experimental  data 
and  showed  a  severe  disagreement  with  the  selectivity  data. 
After  including  the  reverse  water  gas  shift  reaction  and  using 


— 34104\ 


RT 


■ 


Fig.  6.  Comparison  of  experimental  exiting  CO  conversion  and  se¬ 
lectivity  data  with  calculated  numerical  value  from  three  rate  expres¬ 
sions  obtained  this  study.  0.5  g  of  Selectoxo  (Pt-Fe/y-alumina)  catalyst, 
flow  rate  =  167  seem,  FC  =  62-72  mol%,  O2/CO  =  0. 5-2.0,  reaction 
temperature  =  100-300  °C,  1  atm,  no  water  addition. 


the  three  reaction  model,  the  calculated  values  from  the  sim¬ 
ulation  were  almost  exactly  matched  with  the  experimental 
data. 

A  hypothetical  cause  of  decreasing  selectivity  at  higher 
temperatures  might  be  that  the  activation  energy  for  H2  ox¬ 
idation  is  greater  than  the  activation  energy  for  CO  oxida¬ 
tion.  However,  it  was  impossible  to  fit  the  data  accurately  if 
we  constrained  the  activation  energy  in  this  manner.  Rather, 
we  found  the  activation  energy  of  the  H2  oxidation  reaction 
to  be  lower  (one-half)  than  that  for  CO  oxidation.  This  fur¬ 
ther  validates  our  need  to  include  the  WGS  reaction  since 
this  is  the  only  other  reason  why  selectivity  would  decrease 
at  increasing  temperatures. 

4.2.  PROX  reactor  simulations  and  heat  balance  insights 

Using  the  rate  expressions  obtained  above,  a  series  of 
simulations  of  the  PROX  reactor  were  run  based  on  a  1  kW 
fuel  cell/fuel  reformer  system.  Design  parameters  and  op¬ 
erating  conditions  of  the  1  kW  system  are  summarized  in 
the  Table  3.  For  this  system,  11.1  molh-1  of  methanol  must 
be  fed  to  the  methanol  reformer.  If  CO  is  reduced  to  the 
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Table  3 


Design  basis  for  PROX  reactor  and  1  kW  fuel  cell  system 


Design  variables 

Values 

Unit 

PROX  reactor 

CO  inlet  molar  flow  rate 

0.443 

gmol  h- 1 

WGS  reactor 

CO  inlet  molar  flow  rate 

11.1 

gmol  h- 1 

H2  inlet  molar  flow  rate 

22.2 

gmol  h- 1 

CO  target  conditions 

1 

mol% 

Fuel  reformer 

Required  amount  of  hydrogen 

33.3 

gmol  h- 1 

Methanol  feed  rate 

11.1 

gmolh-1 

Fuel  cell  stack 

Target  power 

1 

kW 

Stack  voltage  (50  cells) 

35 

V 

Cell  amperage 

28.57 

A 

Fuel  cell  efficiency,  assumed 

80 

% 

level  of  1  mol%  through  the  WGS  reactor,  then  the  expected 
inlet  concentration  of  CO  to  the  PROX  reactor  would  be 
0.443  gmolh-1. 

The  performance  of  the  PROX  reactor  for  this  basis  is 
shown  in  Fig.  7.  This  Figure  plots  reactor  temperature,  CO 
concentration,  CO  conversion,  and  CO  selectivity  along 
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Fig.  7.  PROX  reactor  simulation  with  various  heat  transfer  cases  using 
50 g  of  Selectoxo  catalyst,  feed  H2  =  33 gmolh-1,  O2/CO  =  1.2,  water 
addition  =  10%,  and  inlet  temperature  =  200  °C,  Ua  =  0,  15,  30,  00 
(Js-1  cm-3  K-1).  Note  that  negative  conversion  and  selectivity  are  not 
shown. 


the  reactor  path  (catalyst  weight)  for  a  feed  temperature  of 
200  °C.  Because  both  CO  oxidation  and  H2  oxidation  are 
exothermic,  the  reactor  temperature  must  be  actively  con¬ 
trolled.  Four  different  heat  exchange  cases  are  considered: 
(1)  perfect  isothermal  operation,  (2)  adiabatic  operation 
and  (3,  4)  actively  cooled  with  two  heat  transfer  rate  co¬ 
efficients.  For  isothermal  operation,  the  CO  conversion 
gradually  reached  95%  at  the  reactor  exit  while  the  CO 
selectivity  remained  almost  constant  at  approximately  40%. 
If  adiabatic  operation  is  assumed,  the  reactor  temperature 
increases  to  a  maximum  372  °C  near  the  middle  of  the 
reactor,  and  then  decreases  slowly  due  to  the  endothermic 
effects  of  the  reverse  WGS  reaction.  For  the  adiabatic  case, 
the  CO  activity  and  selectivity  decrease  rapidly  after  the 
maximum  temperature,  with  a  maximum  CO  conversion  of 
approximately  76%.  The  decreasing  conversion  and  selec¬ 
tivity  are  due  to  the  reverse  WGS  reaction  since  the  oxygen 
has  been  completely  consumed. 

Also  shown  in  Fig.  7  are  two  cases  with  a  sim¬ 
ple  heat  exchange  model  represented  by  Ua  =15  and 
30  Js-1  cm-3  K-1.  The  overall  heat  transfer  coefficient  (U) 
is  combined  with  characteristic  diameter  of  the  reactor  (a 
=  4 ID)  and  it  is  applied  to  the  energy  balance  for  a  steady 
state  non-isothermal  tubular  reactor. 

dT  =  Ua(Tw  -  T)  +  (-AHRx)(-ri) 

dW  J2niCvi  (  j 


-  H2O=0  - H20=7.5% 

■1-120=14%  - 1-120=20% 


Fig.  8.  Exiting  concentration  and  selectivity  with  the  addition  of 
water  versus  reactor  temperature  using  50  g  of  Selectoxo  catalyst, 
H2  =  33gmolh-\  02/CO  =  1.2,  1%  CO. 
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- O2/CO=0.6 - 0.8  . 1.0 

- 1.2  - 1.4 


Fig.  9.  CO  Conversion  and  selectivity  at  various  O2/CO  ratios  ver¬ 
sus  reactor  temperature:  50  g  of  Selectoxo  catalyst,  H2  =  33gmolh_1, 
water  =  10%,  1%  CO. 

For  the  case  of  Ua  =  30  J  s_1  cm-3  K_1,  the  CO  conversion 
and  selectivity  was  higher  than  the  isothermal,  adiabatic  or 
Ua  =  15  J  s-1  cm-3  K-1  cases.  These  results  show  how  im¬ 
portant  it  is  to  properly  design  the  heat  transfer  system  to 
maximize  the  performance  of  the  PROX  reactor. 

4.3.  The  effects  of  water  addition  and  O2/CO  ratio 

Fig.  8  shows  the  effects  of  water  feed  on  CO  conver¬ 
sion  and  selectivity  for  isothermal  operation  at  a  range  of 
temperature.  As  shown  in  the  figure,  the  CO  conversion  is 
enhanced  by  the  addition  of  water  when  the  reactor  tem¬ 
perature  is  above  200  °C.  Also  the  temperature  window 
where  100%  CO  conversion  is  attainable  becomes  wider  as 
more  water  is  added  and  the  selectivity  remains  unchanged 
at  approximately  40%  in  this  window.  This  positive  influ¬ 
ence  of  adding  water  is  consistent  with  our  experimental 
results  and  others.  Manasilp  and  Gulari  [22]  suggested  that 
the  addition  of  water  reduces  the  activation  energy  for  CO 
oxidation.  Also,  they  suggested  three  possibilities  for  water 
enhancement:  (1)  the  water  gas  shift  reaction,  (2)  the  role 
of  the  adsorbed  hydroxyl  group,  or  (3)  the  change  of  Pt 
metal  state.  Our  simulations  show  that  the  reverse  WGS  is 
most  likely  the  cause  of  this  enhancement. 

Fig.  9  shows  the  effects  of  varying  the  O2/CO  feed  ratio. 
To  reach  100%  CO  conversion,  the  O2/CO  ratio  must  be 


greater  than  one.  For  CO  selectivity  higher  O2/CO  ratio 
is  advantageous  at  low  temperatures  while  a  lower  ratio  is 
better  at  higher  temperatures,  due  to  the  competing  effects 
of  the  two  oxidation  reactions. 

In  summary  these  results  suggest  the  best  operation  will 
be  achieved  at  an  O2/CO  ratio  of  1.0-1. 2,  an  isothermal 
reactor  temperature  of  200-220  °C,  and  a  feed  water  content 
above  20%. 


5.  Conclusion 

A  three  reaction  model  consisting  of  CO  oxidation,  H2 
oxidation,  and  the  water  gas  shift  reaction  is  necessary 
to  understand  the  behavior  of  the  PROX  reaction  system. 
Empirical  rate  expressions  for  the  reactions  were  accu¬ 
rately  derived  from  a  large  set  of  experimental  data,  using  a 
commercial  Pt-Fe  catalyst.  Applying  these  empirical  rela¬ 
tionships  in  a  simple  tubular  reactor  model,  it  was  possible 
to  gain  insights  into  the  design  and  operation  of  the  PROX 
reactor.  The  trend  of  decreasing  CO  conversion  and  selec¬ 
tivity  at  higher  temperatures  is  accurately  predicted  to  be 
caused  by  the  reverse  water  gas  shift  reaction  rather  than  a 
difference  in  the  activation  energies  for  CO  oxidation  and 
H2  oxidation.  Also,  it  is  shown  that  adding  water  should  in¬ 
crease  the  performance  of  PROX  reactors.  The  results  of  this 
study  are  expected  to  be  a  critical  part  of  the  overall  design, 
optimization  and  control  of  a  fuel  cell  and  fuel  reformer 
system. 
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